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Application of DNA Barcoding in Taxonomy and Phylogeny: an Individual Case of COI Partial 
Gene Sequencing from Seven Animal Species. Drohvalenko, M., Mykhailenko, A., Rekrotchuk, M., 
Shpak, L., Shuba, V., Trokhymchuk, R., Utevsky, S., Zinenko, O. — A part of the COI mitochondrial 
barcoding gene was sequenced from seven species of different taxonomical groups: Ambystoma 
mexicanum (Amphibia, Ambystomatidae), Darevskia lindholmi, Lacerta agilis exigua (Reptilia, 
Lacertidae), Erinaceus roumanicus (Mammalia, Erinaceidae), Macrobiotus sp. 1 and 2 (Eutardigrada, 
Macrobiotidae) and Cameraria ohridella (Insecta, Gracillariidae). The sequences were compared with 
available sequences from databases and positioned on phylogenetic trees when the taxa had not yet been 
sequenced. The presence of Mexican axolotls in herpetoculture in Ukraine was confirmed. The partial 
COI genes of the Crimean rock lizard and an eastern sub-species of the sand lizard were sequenced. We 
demonstrated the presence of two tardigrade mitochondrial lineages of the Macrobiotus hufelandi group 
in the same sample from the Zeya Natural Reserve in the Far East: one was nearly identical to the Italian 
M. macrocalix, and the other one is similar to M. persimilis and M. vladimiri. We also confirmed the 
presence of the invasive haplotype “A” of the horse chestnut leaf miner in Ukraine, in line with the 
hypothesized route of invasion from Central Europe. 

Key words: cytochrome oxidase subunit one, barcoding, animals, Ukraine, Ambystoma mexicanum, 
Darevskia lindholmi, Lacerta agilis exigua, Erinaceus roumanicus, Macrobiotus sp., Cameraria ohridella. 


Introduction 


The idea of using certain variable parts of a genome as a unique identifier of species has been widely 
recognized since the first publication in 2003 (Hebert et al., 2003). Currently, global initiatives on the barcoding 
of life contain about 6 million records, including 188 000 barcodes obtained from animals (Ratnasingham, 
Hebert, 2007, 2013). With advanced technologies and a drastic drop in the price of sequencing, the wide and 
practical application of DNA barcoding promises to be a simple, fast and precise method to identify species 
from poorly known groups or in case of difficulty with morphological identification. Furthermore, it bridges 
the gap between morphological and molecular identification of taxa, especially in taxonomically diverse groups 
where shortage of qualified experts, known as the taxonomic impediment exists (de Carvalho et al., 2007). In 
spite of great advances of the method, the barcoding of animal biodiversity on Earth is very hard to achieve. 
Areas outside high-income first world countries, remote regions and biodiversity hotspots are poorly covered. 
Hence, many species are still lacking molecular barcodes. We argue that there is a need for more local efforts in 
filling the global list of barcoded species. 
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Ambystoma mexicanum (Shaw, 1789) (Amphibia, Ambystomatidae) also known as the axolotl or the 
Mexican salamander is a neotenic salamander species and is a part of the Ambystoma tigrinum complex. 
As of 2010, wild axolotls in natural habitats in Mexico City were nearly extinct due to urbanization, water 
pollution and introduction of invasive fish species, such as the carp and the tilapia (Zambrano et al., 2010). The 
species is currently listed in CITES (“Appendices I, II and IIT". CITES Secretariat. CITES. Geneva, Switzerland: 
Convention on International Trade in Endangered Species of Wild Fauna and Flora; Checklist of CITES Species 
checklist.cites.org, Ambystoma mexicanum) and is critically endangered in the wild according to the IUCN 
assessment (Zambrano et al., 2010). Since 1998 its density in Lake Xochimilco has decreased drastically from 
thousands per square kilometer to almost zero in 2013 (Voss et al., 2015). At the same time, the axolotl is 
present in large numbers in captivity, as a popular pet and a model organism. Given its endangered status in the 
wild, the captive population may be critical for preserving the species' genetic diversity. Reliable identification 
to avoid confusion of axolotls with tiger salamanders, is also needed. We have sequenced the COI gene from a 
captive individual from Ukrainian herpetoculture in order to confirm its identity. 

Darevskia lindholmi (Szczerbak, 1962) (Squamata, Lacertidae) is a rock lizard species with distribution 
mostly in the Caucasus. Darevskia spp. are famous because ofthe presence of parthenogenetic forms in the genus. 
D. lindholmi was found in Crimea and is the only genus representative in Ukraine. Cytochrome b sequencing 
and protein allozyme analysis (MacCuloch et al., 2000; Tarkhnishvili et al., 2016) show that D. lindholmi is the 
sister species of D. brauneri (Méhely, 1909). The Crimean D. lindholmi has a unique type of Caucasian lizard 
microsatellite allele (CLat IV; along with other types, Ciobanu et al., 2003; Grechko et al., 2006). The source 
and the pattern of colonization of the Crimean Peninsula are not completely clear. Hybridization and multiple 
colonization events, which had led to the origin of this species, were discussed by Kukushkin et al. (2017). 
Obtaining the barcode sequence and clarification of the phylogenetic position of D. lindholmi may help to 
reconstruct this in the future. 

Lacerta agilis exigua Eichwald, 1831 (Squamata, Lacertidae) is one of the 9 subspecies of the sand lizard 
Lacerta agilis Linnaeus, 1758 (Kalyabina et al., 2001; Kalyabina-Hauf and Ananjeva, 2004). L. a. exigua with 
four other subspecies forms Lacerta agilis eastern subspecies group (Kalyabina et al., 2001; Kalyabina-Hauf 
and Ananjeva, 2004; Andres et al., 2014). In eastern and central Ukraine, this subspecies forms a narrow 
contact zone with L. a. chersonensis Andrzejowski, 1832. The contact zone was traced using morphological 
characters (Zinenko et al., 2005). However, no COI gene sequences have been published yet. Mitochondrial 
DNA sequences may be used as additional molecular marker in the hybrid-zone studies. 

The phylogeny of European hedgehog species, relying on the evolution of 125 mitochondrial gene (He 
et al., 2012), depicts Erinaceus europaeus Linnaeus, 1758, Erinaceus concolor Martin, 1837 and Erinaceus 
roumanicus Barrett-Hamilton, 1900 (Mammalia, Erinaceidae) as a monophyletic clade. Despite the fact that 
hedgehogs are common in Europe, the number of barcoding gene sequences available in open databases is 
relatively small. Mitochondrial COI sequences of E. roumanicus and E. concolor have never been published. 
Considering that Ukraine could be a possible contact zone between Eastern European taxa (Bogdanov et al., 
2009) and their distribution is uncertain, rapid and reliable methods of identification could be very useful. 

Water bears or tardigrades (Parachaela, Macrobiotidae) are small, ubiquitous (less than 1 mm) organisms, 
which can tolerate a wide range of abiotic parameters in an inactive stage. The systematics ofthe group is under 
development and a significant number of new species are described every year (Lee et al., 2017; Perez-Pech 
et al., 2017; Kaczmarek & Michalczyk, 2017; Kaczmarek et al., 2017). However, tardigrades are a challenge to 
taxonomists because oftheir small size, unknown magnitude of intraspecific variation and often a lack of experts 
who are able to identify species. Tardigrades, therefore, have a strong need for a unified and reliable method of 
identification and the creation of a good barcode reference library, in which molecular data must be integrated 
with reliable morphological taxonomic knowledge (Cesari et al., 2013). Being small and light in weight, they 
can potentially travel long distances, resulting in their very diffuse, potentially even global distribution (Pilato, 
2001). Recent studies argue that this is not the case, at least at all times (Guil et al., 2009), and that the issue 
could be partly elucidated using information on barcoding genes sequences. 

The horse-chestnut leaf miner Cameraria ohridella Deschka & Dimić, 1986 (Lepidoptera, Gracillariidae) is 
a dangerous pest of horse-chestnut (Aesculus hippocastanum L., Sapindaceae) in Europe. It has been expanding 
the native range in the Balkans since the end of 20th century (Lees et al., 2011) and had colonized Ukraine by 
2007 (Zerova et al., 2007; Guglya, Zinenko, 2008). Among the 25 different mitochondrial haplotypes, only 
three haplotypes took place in the invasion (Valade et al., 2009). Currently, no material from Ukraine has been 
studied, the origin and composition of local lineages of C. ohridella in Ukraine are not known. 

The study was undertaken as an educational project, where students, under faculty guidance, performed 
the necessary manipulations: DNA sampling, depositing or morphological identification of organisms 
followed by the obtaining and submitting DNA barcoding sequences. Thus, we were promoting the concept 
of barcoding life-forms and at the same time facilitating individual research. During the study, we have 
tested the universal primers, optimized reaction conditions and sequenced the barcoding gene COI from 
various (pet, invasive, native) species, which were previously out of the focus of barcoding. Based on the 
obtained sequences, we were able to identify the species or the interspecific lineages using bioinformatics or 
phylogenetic methods. 
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Material and methods 


The samples were conserved in ethanol or taken fresh (vertebrates). Entire larva specimen (e. g. 
C. ohridella) or several specimens (e. g. water bears) were taken. A list of accession numbers of the sequences, 
taxa and sample origin is given in table 1. 

DNA from the tardigrades was extracted in January 2017, one year after sampling. We utilized standard 
methods of tardigrade extraction: a piece of the moss was soaked in water for two hours and the tardigrades 
were picked up from the sample using a stereomicroscope. We also obtained some tardigrade eggs from the 
same sample as they were necessary for species identification. Some of the extracted tardigrades were used 
for DNA analysis. The others from the same sample were used for making permanent slides for subsequent 
microscopy and were mounted in Faure’s medium. 

The preliminary examination of the permanent slides under the microscope had shown the presence of a 
single morphotype of Macrobiotus aff. hufelandi in the text and table (Macrobiotus species belonging to the so 
called *hufelandi group’) tardigrades. We had made several attempts to isolate DNA from individual specimens 
of tardigrades as recommended by Cesari at al., 2011, but none of them was successful. Therefore, we extracted 
DNA from 3 pools of 5 tardigrade specimens torn using a needle. 

Since species of the M. hufelandi group are practically indistinguishable under the magnification of a 
stereomicroscope, our material for the single DNA analysis contained individuals of different species. Later, 
a more thorough examination of the eggs revealed two different species of M. hufelandi group in one sample. 
One of them was later identified as Macrobiotus cf. macrocalix Bertolani & Rebecchi, 1993. The second one was 
Macrobiotus sp. not identifiable to the species level. 

We used high salt DNA extraction method (Aljanabi and Martinez, 1997). Extracted DNA was dissolved 
in water and stored at -18 ° C. Amplification was performed in 25 uL volume, using master mix which con- 
tained DNA, Taq polymerase, dNTPs, primers and 10x buffer. Primers for the COI region were 

LCO1490: GGTCAACAAATCATAAAGATATTG and 

HCO2198: TAAACTTCAGGGTGACCAAAAAAT — shorter by one and two base pairs from 3' end 
than primers which were originally proposed by Folmer et al. (1994). PCR products were obtained using “Ter- 
cik’ and ‘Biometra’ thermocyclers. PCR conditions were optimized by changing annealing temperature from 55 
to 47 °C, until clear single band was visible on 1 % agarose gel electrophoregram. Thermocycling conditions for 
the reaction were as follows: denaturation for 5 min, 36 cycles of denaturation for 60 sec at 94 *C, annealing for 
90 sec at 47 °C and elongation for 60 sec at 72 °C, terminal elongation was set for 5 min at 72 °C. PCR products 
were purified using spin column and then sequenced by commercial sequencing service. 

Chromatograms were manually edited and trimmed. In case of water bears, when several specimens were 
pooled in one sample, the chromatogram contained multiple double peaks, which were first manually coded 
using IUPAC ambiguity codes and then phased using DNASP 5.0 with COI sequences of other Macrobiotus 
spp. retrieved from GenBank. All other samples did not contain any ambiguities and when aligned in a reading 
frame against the COI reference sequences, did not contain any stop codons or indels. 


Table 1. Material and accession numbers of samples sequences 




















: : Date of Accession 
Sample ID Species Locality collection Collector number 
OZ17 Darevskia Ukraine, Crimea, Chufut 2012 Oleg MG458880 
lindholmi Kale Kukushkin 
OZ27 Lacerta agilis Ukraine, Kharkiv Region, 1.07.2016 Oleksandr MG815780 
exigua Zmiiv District, Gaidary Zinenko 
village, Biological station of 
Kharkiv University 
OZ31 Erinaceus Ukraine, Kharkiv Region, 1.07.2016 ‘Tatiana MH329869 
roumanicus Zmiiv District, Gaidary Atemasova 
village, Biological station of 
Kharkiv University 
OZ39 Cameraria Ukraine, Khmelnitskyi 1.08.2016 Oleksandr MH059571 
ohridella Region, Kamiyanets- Zinenko 
Podilskyi 
RF 1.2 Macrobiotus gr. Russia, Amur Oblast, Zeya 2016 Anna Suvorova MH042875 
hufelandi (?) Nature Reserve 
ARI Ambystoma Ukraine, Kharkiv, captive 2017 Roman MH791446 


mexicanum bred Trokhymchuk 
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Table 2. Parameters of phylogenetic tree reconstruction 


























Best substitution ; 
Taxa names model under BIC LnL | MrBayes settings 
Lacerta agilis HKY -2392.312 lset applyto = (all) nst = 6 rates = invgamma; mcmc 
exigua and ngen = 1000000 
Darevskia 
lindholmi 
Macrobiotus GTR+G+I -3870.424 lset applyto = (all) nst = 6 rates = invgamma; mcmc 
hufelandi ngen = 2000000 samplefreq = 100 printfreq = 100 
group diagnfreq = 1000 
Erinaceus Pos 1 and 2 — -741.046 charset Subset] = 1-286\3 2-286\3; charset Subset2 = 
roumanicus K80+I; 3-286\3; partition PartitionFinder = 2:Subset1, Subset2; 
pos 3 — HKY+G set partition = PartitionFinder; lset applyto = (1) 


nst = 2 rates = propinv; prset applyto = (1) 

statefreqpr = fixed(equal); Iset applyto = (2) nst = 2 
rates = gamma; prset applyto = (all) ratepr = variable; 
unlink statefreq = (all) revmat = (all) shape = (all) 
pinvar = (all) tratio = (all); unlink brlens = (all); meme 
ngen = 1000000; 





Sequences were exported in the FASTA format and similar sequences were retrieved from GenBank 
and BOLD databases (Ratnasingham and Hebert, 2007). To resolve phylogenetic relationships within studied 
genera, COI sequences from every species within the group were taken, with the addition of one species outside 
the genus as an outgroup. 

The sequences were initially identified using BLAST search. If BLAST search was able to find identical or 
closely identical sequences (in the case of Ambystoma, Cameraria samples), further analysis was not performed. 
In the case of Erinaceus sp., Macrobiotus sp., Lacerta sp. and Darevskia sp., samples were identified at the 
genus level using BLAST and reference sequences from each genus were included into further analyses. The 
reference and query sequences of the genus were aligned using ClustalW in BioEdit 7.1.9 (Hall, 1999) and Mega 
7.0.26 (Kumar et al., 2016). Optimal substitution model was calculated in Mega 7.0.26 and Partition Finder 
2.0. Topology of phylogenetic trees was constructed in MrBayes 3.2.6 (Bayesian inference - BI) or Mega 7.0.26 
(Maximum Likelihood, ML). Trees topology was shown in FigTree 1.4.2. Parameters of tree construction are 
given in table 2. 


Results 


The sequence (337 bp) of the captive bred axolotl did not show any differences from 
the published GenBank sequences AY65999 land AJ584639 of the axolotl. It had two mu- 
tations, which differentiate it from Ambystoma andersoni Krebs and Brandon, 1984. 

According to our 539 base pairs sequence of COI gene, Darevskia lindholmi is in a 
monophyletic clade with other available sequences of lizards from the Darevskia genus 
(D. unisexualis and D. valentini), Eremias and Ophisops (fig. 1). Podarcis, Lacerta and Apa- 
thya genera are shown to be sister clades. The partial COI (630 bp) sequence from a mor- 
phologically identified L. a. exigua was published for the first time. Although at least two 
large scale phylogeographic studies were published (Kalyabina-Hauf and Ananjeva, 2004; 
Andres et al., 2014), the authors used different mitochondrial loci. Ukrainian COI sequence 
was clustered with other representative of the species L. agilis (fig. 2), but not identical to 
them, as it would be expected since GenBank sequences belong to the Central European 
subspecies group L. a. agilis, in contrast to L. a. exigua from the Eastern subspecies group 
(fig. 1). 

The short E. roumanicus COI fragment of 288 bp from Kharkiv Region is clustered first 
with additional sequences of this species and then with sister species E. europeus. E. rou- 
manicus and E. europeus are nested within a larger cluster which includes E. amurensis 
sequences (posterior probability 1) (fig. 2). 

After phasing the Tardigrada sequences we have obtained two 480 base pairs sequenc- 
es of Macrobiotus sp. There were 87 ambiguities in the complete sequence submitted to the 
GenBank, and 80 ambiguities in the final sequence used for further analyses. After phasing, 
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KP993006 


KP993002 Erinaceus amurensis 
KP993001 





KP993000 
xeess Erinaceus europeus 
MH329869 

ky7sas0s | Erinaceus roumanicus 
SKMZM654 

KY605342 

KY605341 

«soso | Neotetracus sinensis 
KY605339 

HM031766 


Hmo31765 | Neohylomys hainanensis 


HM031764 


0.87 


0.009 
Fig. 2. Phylogenetic position of E. roumanicus, Bayesian inference phylogenetic tree. Sequences obtained by us 
are written in bold. 


all ambiguities were resolved with 1.0 probability assignment. The two resulting sequenc- 
es cluster either with M. macrocalix FJ176217 from Italy (p-distance to the reference se- 
quence = 0.009) or with M. vladimiri JX 683811 from Portugal and M. persimilis EU244608 
from Germany (p-distance = 0.020). 

The sequence of C. ohridella from Kamyanets-Podilsky contained 505 base pairs and 
had one substitution comparing to invasive haplotype A, found in the Balkans, Central and 
Western Europe (Valade et al., 2009; Lees et al., 2011). 


Discussion 


Short fragment of the COI gene has confirmed the identity of the captive bred of A. 
mexicanum from Ukraine. Recent radiation of the Ambystoma genus, small length of the 
sequence, which contains only a few substitutions and is not species-specific, lead to the 
conclusion that one mitochondrial gene is not sufficient to trace the origin or define any 
intraspecific relations, but can help differentiate between the most closely related sister 
species A. mexicanum and A. andersoni. 

We also have shown that the universal primers and COI sequences from E. roumanicus 
and L. a. exigua could be used for their rapid identification. Furthermore, it can be used to 
identify the maternal lineage in the individuals from the hybrid zones between L. a. exigua 
and L. a. chersonensis (Zinenko et al., 2005) or between Erinaceus species in Eastern Europe 
and the Caucasus. In both cases, sampling localities are within the distribution ranges of 
particular taxa and could be unambiguously associated with L. a. exigua (see Kalyabina- 
Haufand Ananjeva, 2004; Zinenko et al., 2005) or E. roumanicus (He et al., 2012), therefore, 
they can be used as a reference. 
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aa Macrobiotus persimilis Germany EU244608 
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m Macrobiotus hufelandi isolate KC146124 
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Fig. 3. Phylogenetic position of Macrobiotus sp., Bayesian inference phylogenetic tree. Sequences obtained by 
us are written in bold. 
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The Lacertidae tree with the COI sequence of D. lindholmi was compared to a previously 
described Lacertidae phylogeny (Kapli et al., 2011). The latter included the representatives 
of 40 genera and the sequences of 6 genes (including COI). Our tree differs from it by 
topology on the level of genera. Statistical support of nodes in the tree provided by this study 
was lower due to too short sequence of the COI fragment, which lacked phylogenetically 
informative substitutions. Our tree reconstruction had high support at the level of genera, 
but higher systematics levels above the genus level appeared to be more tangled compared to 
other authors. Hence, we should use more loci and longer sequences to investigate relations 
among Darevskia as well as other genera of lacertid lizards. Intraspecific relations could be 
revealed using denser sampling of closely related species from the Caucasus, which were 
not available at the time of the study. 

Our chosen sampling strategy and sample processing of Tardigrada was successful, 
however we would recommend this approach only for those localities with well-known fauna 
of Tardigrada or in metabarcoding studies. Our sampling area was probably too remote to 
host the exact lineages of European Macrobiotus species, whose sequences from GenBank 
were taken for comparison. Nevertheless, the genetic distance between the obtained phased 
sequences and the nearest sequences on the tree were rather small, confirming the high 
dispersal abilities of water bears. Finally, the presence of two mitochondrial lineages from 
different species ofthe same genus in one environmental sample is rather interesting. On this 
stage we are too far away from understanding the ecology of these species. Mitochondrial 
markers alone are not enough to understand ifthey truly coexist as separate taxa or whether 
gene flow is possible between them (fig. 3). 

The finding of a widespread invasive haplotype of C. ohridella in Ukraine confirms 
the high invasive potential of this particular lineage, which has been shown to substitute 
even other non-invasive haplotypes in the natural range of the species in the Balkans (Lees 
et al., 2011). Invasion of this haplotype to Ukraine is consistent with the vehicle-mediated 
dispersal of C. ohridella and the Central European origin of the species in Ukraine (Zerova 
et al., 2007). 

Standard methods of DNA barcoding could be used to answer various research 
questions. We have successfully utilized them to elucidate taxonomic or intraspecific 
position of A. mexicanum, Macrobiotus sp. and C. ohridella. Newly obtained partial gene 
COI sequences from L. a. exigua, D. lindholmi and E. roumanicus could be used to study 
their distribution patterns, hybridization and phylogeography in Eastern Europe. 


We would like to thank Oleksii Solodyankin for his continuous help in the laboratory establishment — 
Meike Konderman for providing protocols, Viktoriia Komarista for further optimization, Andrii Utevsky, 
Ulrich Joger and Miguel Vences for donating equipment, Anna Suvorova for sampling and Yevgen Kiosya for 
morphological identification of water bears, Oleg Kukushkin and Tatiana Atemasova for samples collection, 
numerous students of Biology faculty of Karazin University for working on optimization of the protocols 
and Oleksii Yanchukov for the valuable help with editing manuscript. This study was supported by charitable 
organization Ukrainian Fulbright Circles small grant, holded by Oleksandr Zinenko in 2016 and partly 
financed within the research project by Serge Utevsky, # 011740048360 “Biological challenges and threats due 
to migrations and invasions" funded by Ukrainian ministry of science and education. Mykola Drohvalenko, 
Anastasiia Mykhailenko, Mariia Rekrotchuk, Liubov Shpak, Vladislav Shuba and Roman Trokhymchuk 
contributed to the paper equally and are listed in alphabetical order. 


References 


Aljanabi, S. M., Martinez, I. 1997. Universal and rapid salt-extraction of high quality genomic DNA for PCR- 
based techniques. Nucleic Acids Res, 25, 4692-4693. 

Andres, C., Franke, F, Bleidorn, C, Bernhard, D, Schlegel, M. 2014. Phylogenetic analysis of the Lacerta agilis 
subspecies complex. Syst Biodivers, 12, 43-54. 


Application of DNA Barcoding in Taxonomy and Phylogeny... 383 


Bogdanov, A. S., Bannikova, A. A., Pirusskii, Y. M., Formozov, N. A. 2009. The First Genetic Evidence of 
Hybridization between West European and Northern White-breasted Hedgehogs (Erinaceus europaeus 
and E. roumanicus) in Moscow Region. Biol Bull Russ Acad Sci, 36 (6), 647—651. 

Cesari, M, Guidetti, R, Rebecchi, L, Giovannini, I, Bertolani, R. 2013. A DNA barcoding approach in the study 
of tardigrades. J Limnol, 72, 23. 

Ciobanu, D. G., Grechko, V. V., Kramerov, D. A., Darevsky, I. S. 2003. A New Subfamily of the Satellite DNA, 
CLsatIV, of the Lizard Darevskia lindholmi (Sauria, Laceridae): Structure and Evolution. Dokl Biochem 
Biophys, 392 (1-6), 263-267. 

De Carvalho, M. R., Bockmann, F. A., Amorim, D. S., Brandao, C. R. F., de Vivo, M., de Figueiredo, J. L., 
Britski, H. A., de Pinna, M. C. C., Menezes, N. A., Marques, F. P. L., Papavero, N., Cancello, E. M., 
Crisci, J. V., McEachran, J. D., Schelly, R. C., Lundberg, J. G., Gill, A. C., Britz, R., Wheeler, Q. D., 
Stiassny, M. L. J., Parenti, L. R., Page, L. M., Wheeler, W. C., Faivovich, J., Vari, R. P., Grande, L., 
Humphries, C. J., DeSalle, R., Ebach, M. C., Nelson, G. J. 2007. Taxonomic Impediment or Impediment 
to Taxonomy? A Commentary on Systematics and the Cybertaxonomic-Automation Paradigm. Evol Biol, 
34, 140-143. 

Folmer, O., Black, M., Hoeh, W., Lutz, R., Vrijenhoek, R. 1994. DNA primers for amplification of mitochondrial 
cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotech, 3 (5), 294-299. 

Grechko, V. V., Ciobanu, D. G., Darevsky, I. S., Kosushkin, S. A., Kramerov, D. A. 2006. Molecular evolution 
of satellite DNA repeats and speciation of lizards of the genus Darevskia (Sauria: Lacertidae). Genome, 49 
(10), 1297-1307. 

Guil, N., Sanchez-Moreno, S., Machordom, A. 2009. Local biodiversity patterns in micrometazoans: Are 
tardigrades everywhere? Syst Biodivers, 7 (3), 259-268. 

Guglya, Yu., Zinenko, O. 2008. New Data on Expansion of the Leaf-Mining Moth, Cameraria ochridella 
(Lepidoptera, Gracillariidae), in Ukraine. Vestnik Zooogiil, 42 (3), 220. 

Hall, T. A. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for 
Windows 95/98/NT. Nucl. Acids. Symp. Ser., 41, 95-98. 

He, K., Chen, J. H., Gould, G. C., Yamaguchi, N., Ai, H-S., Wang, Y-X., Zhang, Y-P., Jiang, X-L. 2012. An 
estimation of Erinaceidae phylogeny: a combined analysis approach. PLoS One 7: e39304. 

Hebert, P. D. N., Cywinska, A., Ball, S. L., de Waard, J. R. 2003. Biological identifications through DNA 
barcodes. Proc Biol Sci, 270 (1512), 313-321. 

Kaczmarek, L., Gawlak, M., Bartels, P. J., Nelson, D. R., Roszkowska, M. 2017. Revision of the Genus 
Paramacrobiotus Guidetti et al., 2009 with the Description of a New Species, Re-Descriptions and a Key. 
Ann Zool, 67, 627-656. 

Kaczmarek, L., Michalczyk, L. 2017. The Macrobiotus hufelandi group (Tardigrada) revisited. Zootaxa, 4363 
(1), 101-123. 

Kalyabina, S. A., Milto, K. D., Ananjeva, N. B., Legal, L., Joger, U., Wink, M. 2001. Phylogeography and 
systematics of Lacerta agilis based on mitochondrial cytochrome b gene sequences: first results. Russ J 
Herpetol, 8, 149-158. 

Kalyabina-Hauf, S. A., Ananjeva, N. B. 2004. Phylogeography and intraspecies structure of wide distributed 
sand lizard, Lacerta agilis L. 1758 (Lacertidae, Sauria, Reptilia). Trudy Zoologicheskogo Instituta, 302, 
1-105. 

Kapli, P., Poulakakis, N., Lymberakis, P., Mylonas, M. 2011. A re-analysis of the molecular phylogeny of 
Lacertidae with currently available data. Basic Appl Herpetol, 25, 97-104. 

Kukushkin, O. V., Doronin, I. V., Tuniyev, B. S., Ananjeva, N. B., Doronina, M. A. 2017. Introduction of 
Amphibians and Reptiles at the Caucasus and the Crimea: an Overview and Some Actual Data. Curr Stud 
Herpetol, 17, 157-197. 

Kumar, S., Stecher, G., Tamura, K. 2016. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for 
bigger datasets. Mol Biol and Evol, 33 (7), 1870-1874. 

Lee, J., Rho, H. S., Chang, C. Y. 2017. Taxonomic Study of Marine Tardigrades from Korea III. A New Species 
of the Genus Orzeliscus (Heterotardigrada, Halechiniscidae). Animal Syst Evol Divers, 33, 26-32. 

Lees, D. C., Lack, H. W., Rougerie, R., Hernandez-Lopez, A. Raus, T., Avtzis, N., Augustin, S., Lopez- 
Vaamonde, C. 2011. Tracking origins of invasive herbivores through herbaria and archival DNA: the case 
of the horse-chestnut leaf miner. Front Ecol Environ, 9, 322-328. 

MacCuloch, R. D., Fu, J., Darevsky, I. S., Murphy, R. W. 2000. Genetic evidence for species status of some 
Caucasian rock lizards in the Darevskia saxicola group. Amphibia-Reptilia, 21, 169-176. 

Pérez-Pech, W. A., Anguas-Escalante, A., Cutz-Pool, L. Q., Guidetti, R. 2017. Doryphoribius chetumalensis sp. 
nov. (Eutardigrada: Isohypsibiidae) a new tardigrade species discovered in an unusual habitat of urban 
areas of Mexico. Zootaxa, 4344 (2), 345-356. 

Pilato, G., Binda, M. G. 2001. Biogeography and limno-terrestrial tardigrades: are they truly incompatible 
binomials? Zool Anz, 240, 511-516. 

Ratnasingham, S., Hebert, P. D. N. 2007. BOLD: The Barcode of Life Data System (www.barcodinglife.org). Mol 
Ecol Notes, 7, 355-364. 


384 M. Drohvalenko, A. Mykhailenko, M. Rekrotchuk, L. Shpak, V. Shuba, R. Trokhymchuk et al. 


Ratnasingham, S., Hebert, P. D. N. 2013. A DNA-Based Registry for All Animal Species: The Barcode Index 
Number (BIN) System. PLoS ONE 8, e66213 

Tarkhnishvili, D., Gabelaia, M., Mumladze, L., Murtskhvaladze, M. 2016. Mitochondrial phylogeny of the 
Darevskia saxicola complex: two highly deviant evolutionary lineages from the easternmost part of the 
range. Herpetol J, 26, 175-182. 

Valade, R., Kenis, M., Hernandez-Lopez, A, Augustin, S., Mari Mena, N., Magnoux, E., Rougerie, R., Lakatos, F., 
Roques, A., Lopez-Vaamonde, C. 2009. Mitochondrial and microsatellite DNA markers reveal a Balkan 
origin for the highly invasive horse-chestnut leaf miner Cameraria ohridella (Lepidoptera, Gracillariidae). 
Mol Ecol, 18, 3458-3470. 

Voss, S. R., Woodcock, M. R., Zambrano, L. 2015. A tale of two axolotls. BioScience, 65, 1134-1140. 

Zambrano, L., Valiente, E., Vander Zanden, M. J. 2010. Food web overlap among native axolotl (Ambystoma 
mexicanum) and two exotic fishes: carp (Cyprinus carpio) and tilapia (Oreochromis niloticus) in 
Xochimilco, Mexico City. Biol Invasions, 12, 3061-3069. 

Zambrano, L., Reidl, P. M., McKay, J., Griffiths, R., Shaffer, B., Flores-Villela, O., ParraOlea, G., Wake, D. 2010. 
Ambystoma mexicanum. In: IUCN 2010. The IUCN Red List of Threatened Species. Version 2010. DOI: 
10.2305/IUCN.UK.2010-2.RLTS.T1095A3229615.en. 

Zerova, M. D., Nikitenko, G. N., Narolsky, N. B., Gershenson, Z. S., Sviridov, S. V., Lukash, O. V., 
Babidoritsh, M. M. 2007. Horse-chestnut Leaf miner, Cameraria ohridella, in Ukraine. Kiev, Ukraine, 
Schmalhausen Institute of Zoology [In Russian]. 

Zinenko, O. L, Drabkin, P. L., Rudyk, O. M. 2005. Contact zone between two subspecies of the Sand lizard: 
Lacerta agilis exigua Eichw., 1831 and Lacerta agilis chersonensis Andr., 1832 in three regions of the Left- 
bank Ukraine. In: Ananjeva, N., Tsinenko, O., eds. Herpetologica Petropolitana. Proceedings of the 12th 
Ordinary General Meeting Societas Europea Herpetologica, August 12-16, 2003, St. Petersburg, Russian 
Federation: SEH, 109-112. 


Received 30 January 2019 
Accepted 5 March 2019 


